Acceleration disturbances and requirements for 
ASTROD I 



Sachie Shiomi 1 and Wei-Tou Ni 1 ' 2,3 

department of Physics, National Tsing Hua University, Hsinchu, Taiwan 30013 
ROC 

2 Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210008 
China 

3 National Astronomical Observatories, Chinese Academy of Sciences, Beijing 
100012 China 



Abstract. 

ASTRODynamical Space Test of Relativity using Optical Devices I (ASTROD 
I) mainly aims at testing relativistic gravity and measuring the solar-system 
parameters with high precision, by carrying out laser ranging between a spacecraft 
in a solar orbit and ground stations. In order to achieve these goals, the magnitude 
of the total acceleration disturbance of the proof mass has to be less than 
10 -13 ms -2 Hz -1 / 2 at O.lmHz. In this paper, we give a preliminary overview 
of the sources and magnitude of acceleration disturbances that could arise in the 
ASTROD I proof mass. Based on the estimates of the acceleration disturbances 
and by assuming a simple control-loop model, we infer requirements for ASTROD 
I. Our estimates show that most of the requirements for ASTROD I can be relaxed 
in comparison with Laser Interferometer Space Antenna (LISA). 



PACS numbers: 04.80.Nn, 04.80.-y, 95.55.Ym 

1. Introduction 

ASTRODynamical Space Test of Relativity using Optical Devices (ASTROD) P E] 
aims at testing relativistic gravity, measuring the solar-system parameters with high 
precision and detecting gravitational waves from massive black holes and galactic 
binary stars. The concept of ASTROD is to put two spacecraft in separate solar 
orbits and carry out laser interferometic ranging with Earth reference stations (e.g. a 
spacecraft at the Earth-Sun L1/L2 points). A simple version of ASTROD, ASTROD 
It, has been studied as the first step to ASTROD E]. ASTROD I employs 
one spacecraft in a solar orbit and carries out interferometric ranging and pulse 
ranging with ground stations. The main scientific goals of ASTROD I are to 
test relativistic gravity and the fundamental laws of spacetime with three-order-of- 
magnitude improvement in sensitivity and to improve the solar, planetary and asteroid 
parameter determination by 1 to 3 orders of magnitude. The technological goal of 
ASTROD I is to prepare for the ASTROD mission. 

J ASTROD I was previously referred to as Mini- ASTROD (for instance in [3*|'l. 
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The acceleration disturbance goal of the ASTROD I proof mass is 
10 -13 ms~ 2 Hz -1 / 2 at frequency v of O.lmHz. The power spectral density of the al- 
lowed level of the acceleration noise is shown in figure^ Assuming a 10 ps one-way tim- 
ing accuracy (3 mm ranging accuracy) and the acceleration noise of 10 -13 m s -2 Hz -1 / 2 
at frequency of about O.lmHz, a simulation for 400 days (350-750 days after launch) 
showed that ASTROD I could determine the relativistic parameters 7 and /3, and the 
solar quadrupole parameter J2 to levels of 10~ 7 , 10~ 7 and 10 -8 , respectively (SJ. In 
the simulation, (i) the timing noise is modeled as Gaussian random noise; (ii) un- 
known acceleration noise is modeled to have Gaussian random magnitude with zero 
mean and with standard deviation 10 -15 ms~ 2 and to have its direction changed ran- 
domly every 4h (equivalent to 10~ 13 ms~ 2 Hz -1 / 2 for v ~ O.lmHz assumed as the 
requirement of the drag-free system) and (iii) five range points are taken each day (at 
0.2 d interval). Longer term systematic effects will be studied in a future paper. This 
simulation agrees with the scientific goals of ASTROD I. The timing uncertainty of 
event timer reaches 3ps in satellite laser ranging at present. Space qualified versions 
of similar accuracy are under development. For a ranging uncertainty of 3 mm in a 
distance of 3 x 10 11 m (2 AU), the laser/clock frequency needs to be known to one part 
111 10 14 . This can be set as a requirement of the space laser/clock or a requirement for 
laser frequency monitoring through ground clock and modeling. As to ground station 
jitter, monitoring to an accuracy of 3 mm is required and can be achieved. The atmo- 
spheric effects on laser propagation will be monitored and subtracted to mm-level by 
using 2-color (2-wavelength) ranging (one color for pulse ranging and one for interfer- 
ometric ranging). These measurement uncertainties are not cumulative in the range 
determination while the acceleration disturbances accumulate in time in the geodesic 
deviations. In order to achieve the acceleration disturbance goal, a drag-free control 
system using capacitive sensors will be employed. 

In order to decide on detailed designs of the accelerometer for ASTROD I, we 
have to know the sources and magnitude of acceleration disturbances that could arise 
in the accelerometer. In this paper, we carry out analyses, mainly based on existing 
literature on acceleration disturbances for other gravitational missions (LISA J7| |Hj 
and the LISA Pathfinder El)) to give a preliminary overview of the acceleration 
disturbance to the ASTROD I proof mass. Based on the analyses, we infer some of 
the requirements for the designs of the ASTROD I payload and spacecraft. Also, we 
compare parameter values we have assumed in the analyses with those for LISA to 
confirm the feasibility of the requirements for ASTROD I. 

First, we will give an overview of the ASTROD I configuration (section^ and 
the control-loop model we assumed (section |3J). Then, we will estimate the magnitude 
of acceleration disturbances and requirements for ASTROD I in sections 0] to [HI and 
compare the requirements of ASTROD I with LISA in section OH 

2. ASTROD I spacecraft configuration 

The ASTROD I spacecraft has a cylindrical shape with diameter 2.5 m and height 2 m. 
Its cylindrical side is covered by solar panels. The cylindrical axis is perpendicular 
to the orbit plane and the telescope is set to point toward a ground laser station. 
The effective area of receiving sunlight is about 5 m 2 and it can generate power that 
is larger than 500 W. The total mass of the spacecraft is about 350 kg and that of 
payload is 100-120 kg (see f° r more detailed descriptions of the configuration). 

The orbit distance from the Sun varies from about 0.5 AU to 1 AU (figure 2 of fTU]l. 
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Figure 1. Acceleration noise spectral density requirements for ASTROD I. The 
LISA Technology Package (LTP) requirements |§] and LISA requirements [2] are 
illustrated in the figure for comparison. 

The proof mass (m p — 1.75 kg) is a rectangular parallelepiped (50 x 50 x 
35mm 3 )§ made from Au-Pt alloy (density p = 2 x 10 4 kgm~ 3 ). The six sides of 
the proof mass are surrounded by electrodes mounted on the housing for capacitive 
sensing and control. The gap between each side of the proof mass and the opposing 
electrode is 2 mm. Assumed values for the capacitance and voltages are listed in table 

m 

3. Control-loop model 

Various acceleration disturbances would act on the proof mass in different ways. In 
order to infer how these different kinds of acceleration disturbances would contribute to 
the total acceleration disturbance of the proof mass, we tentatively assume a simple 
single-mass and single-axis control- loop model | 1 II I12| . The diagram is shown in 
figure The relative difference between the displacement disturbance amplitudes of 
the proof mass and of the spacecraft, X ps = X p — X s , with position readout noise X nr , 
is measured by the position displacement sensor. The output of this sensor is converted 
to acceleration disturbance f r , by a transfer function R (= w^). This acceleration is 
supplied to the thruster and the output acceleration disturbance with thruster noise N t 
is applied to the spacecraft. The spacecraft also experiences acceleration disturbances 
by coupling to the proof mass with a coupling constant K and external environmental 

§ This is the current design of the proof mass. A cylindrical shape is also considered as an alternative 
design of the proof mass. 
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Figure 2. Control loop model. SC and PM denote spacecraft and proof mass, 
respectively. 



disturbances (/ ns ) that work directly on the spacecraft. The total acceleration of the 
spacecraft f s is converted to the position noise X s with a transfer function S. 

We consider only the sensitive axis of the proof mass without actuation. The 
proof mass would experience disturbances by spacecraft-proof mass coupling and 
environmental disturbances f np (see section |5J. The total acceleration disturbance 
of the proof mass f p is converted to the displacement disturbance X p with a transfer 
function P. 

From this control-loop model, we obtain the following linear loop equations 

Q31G3]: 



fr — R(X ps + X nr ) 


(1) 


ft =T(f r + N t ) 


(2) 


fs = fns + ft~ (m p /M sc )f ps 


(3) 


fp — fnp "i" fps 


(4) 


fps KXp S 


(5) 


X s = Sf s 


(6) 


Xp — Pfp 


(7) 



By solving these equations for X p , and assuming S = P = u> 2 and u = STR (= 
TujftUJ^ 2 ; T is nominally 1 and its effects are absorbed in R and Nt), we obtain (e.g. 
[Hid Eld): 



fp w X nr (-K) + f np + (f ns + TN t )Kui~ 



(8) 
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where w — lixv. This acceleration disturbance has to be less than the acceleration 
noise goal of 10 -13 ms -2 Hz -1 / 2 at v = O.lmHz. 

We will estimate the values of the direct acceleration disturbances of the 
spacecraft (f ns and the thruster noise TN t ) in section 0] and of the proof mass (/„ p ) 
in sections and El and the stiffness K in section We will discuss the requirements 
for X nr and u in section |H1 

4. Direct acceleration disturbances of the spacecraft 

The spacecraft would be affected by environmental disturbances that stem from, for 
example, solar radiation pressure, solar wind and micrometeorite impacts. Among 
these sources of disturbances, solar radiation pressure is considered to be the major 
contributor to the acceleration disturbances (section 7 of 0). The contribution from 
solar wind might be comparable to radiation pressure, but the spectral behavior of 
the solar wind is not well known. 

By assuming a perfectly reflecting surface of the spacecraft, acceleration noise 
caused by fluctuation in solar irradiance SWo is fns.srp = | (A sc 5Wo/M sc c) , where 
A sc is the area of the spacecraft facing the Sun, M sc is the mass of the spacecraft 
and c is the speed of light in vacuum. From the data of the VIRGO experiment on 
SOHO [To), fractional fluctuation in solar irradiance is (Wo/Wo 2.8 x 10~ 3 Hz~ 1//2 
at O.lmHz (figure 6 of US]). Assuming that SWq/Wq is the same at 0.5 AU, and 
taking a total irradiance of 5500 Wm -2 , the area of 5 m 2 and a mass of 350 kg, we 
obtain /„ s , srp = 9.8 x 10~ 10 ms" 2 Hz^ 1 / 2 . 

The impact rate of 1-ng meteorites on the ASTROD I spacecraft (its surface 
area is 25.5m 2 ) is about 7 events per day or about 0.08mHz||. The average velocity 
of meteorites is 18kms _1 in an Earth frame near the Earth (figure 7.2-30 of 8 ). 
An impact of 1-ng meteorite with the average velocity of 18kms _1 on the surface 
of the spacecraft (350 kg) without reflection produces a linear velocity increment of 
about 5 x 10 _11 ms _1 . Smaller meteorites have larger flux (figure 7.2-29 of 8), 
but their impacts on velocity changes are smaller. Therefore, the contribution to the 
acceleration disturbances at 0.1 mHz seems insignificant, in comparison with the effect 
of solar radiation pressure. However, because such discrete changes could directly 
affect the ASTROD I experiment, the impact effects have to be studied carefully in 
detail. 

In addition to the acceleration noise from the environmental disturbances, 
the spacecraft would suffer from the thruster noise. A force fluctuation of 
10/^NHz -1 / 2 in thruster corresponds to the acceleration disturbance of TNt — 2.8 
x logins" 2 Hz" 1 / 2 . 

Therefore, the rss (root-sum-square) of these acceleration disturbances is: 

f ns + TN t « 2.8 x 10~ 8 m s~ 2 Hz~ 1/2 (9) 

at O.lmHz, dominated by the thruster noise. As per 0, LISA requires a force 
fluctuation to be less than a few /iNHz -1 / 2 . However, recent studies for the LISA 
Pathfinder indicate that force noise up to O.lmNHz -1 / 2 at O.lmHz can be tolerated 
by increasing the gain JSj ; the LISA study team is designing high gain loops and the 
requirement has been relaxed^. 

|| This rate was tentatively estimated from the meteorite flux/mass distribution used by the LISA 

study team (figure 7.2-29 of 0). 

1 This is a comment from one of the referees. 
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5. Direct acceleration disturbances of the proof mass 

Direct proof-mass acceleration disturbances can be classified into two categories 
depending on their origins environmental disturbances (faep) and proof-mass 

sensor back-action acceleration disturbances (fabp)- The former includes disturbances 
related to magnetic effects (f m i, f m 2, fm3, fh\ and Jl2), impact effects (by cosmic 
ra y (/ c ) and residual gas (f r g)), temperature dependent effects (radiometric and 
outgassing effects (f re and f og ) and thermal radiation pressure (fa)) and gravity 
gradients caused by thermal distortion of the spacecraft (f gg ). The latter is originated 
from voltage fluctuations (fa and fa), charge fluctuations (fa and fa), readout 
electronics (fa), patch field voltages (f pe ) and thermal voltage noise by dielectric 
losses (f M ). 

Parameter values and physical constants used for the estimations are listed in 
tables n and respectively. Table provides a summary of the expressions used 
to estimate the direct proof-mass acceleration disturbances and the estimated values. 
The disturbances noted as fat, f m2 , Jli, fhi, fa frg, fa, fa, f gg and fa-fa in 
table correspond to Ai-Aq and Ag,-Au in ^3], respectively" 1 ". We will briefly explain 
each disturbance below. 

5.1. Magnetostatic interaction 

The lowest order force on a proof mass with a magnetic moment M p in an external 
magnetic field B is given by F m = V(M P ■ B) 17 . The magnetic moment of the 
proof mass is a vector sum of the remanent moment M r and the induced moment |18| : 
M p = M r + (x m Vp/f-o)B, where Xm and V p are magnetic susceptibility and the volume 
of the proof mass, respectively; fio = 1.26 x 1CP 6 N A~ 2 is the permeability of vacuum. 
The external magnetic field would be given by the superposition of the interplanetary 
magnetic field Bi p and a local magnetic field B sc . Dominant terms in acceleration of 
the proof mass due to the induced magnetic moment produce acceleration disturbances 
noted fmi and f m i in table 01 where £ rn is a scaling factor for possible suppression by 
magnetic shielding. 

The average interplanetary magnetic field at 1 AU from the Sun varies from 1CP 9 
to 3.7 x 1(T 8 T dnj. Ulysses data obtained near 1 AU from the Sun (figure 9 of 
20 ) showed a v~ 2 / 3 dependence of the variation in the interplanetary magnetic field 
and 5B ip can be inferred to be about 10~ 7 THz -1 / 2 at O.lmHz. As the behavior 
of 6 B ip at 0.5 AU from the Sun is uncertain, we use a somewhat higher value of 
SBi p = 4 x 10 _7 THz~ 1 / 2 (see footnote in section I5~2*|) . According to the formulation 
studies for LISA and the implementation work of the LISA Pathfinder, batteries and 
micro-thrusters are the principal suspects of the origins of the local magnetic field in 
LISA* . We need elaborate modeling works to estimate the magnitude of local magnetic 
field. Here, we use the same values used in analyses for LISA ^]|^ : B sc w 8 x 10~ 7 T, 
SB SC = lO^THz- 1 / 2 and \VB SC \ w 3 B sc r^ « 3 x lO^Tm" 1 (where r m = 0.75m). 
Silvestri et al. have reported that the magnetic susceptibility of five samples ranged 
from —2.8 x 10 -5 to —2.1 x 10~ 5 for two of them without traceable iron contamination 
and from +1.1 x 10~ 5 to +8.8 x 10~ 5 for the rest samples with a trace of iron |22| . 

+ The acceleration disturbance due to laser photon radiation pressure, noted as Ar in 1141. would not 
arise in ASTROD I, as the laser beam is not to be injected directly on the surface of the ASTROD I 
proof mass. 

* This is a comment from one of the referees. 
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Table 1. Parameter values used in the acceleration noise estimates 

v = 0.1 mHz ASTROD I 



Proof Mass (PM) 

Mass: m p [kg] 1.75 

Density: p [kg m~ 3 ] 2 X 10 4 

Cross section: A p [m 2 ] 0.050 X 0.035 

Temperature: T p [K] 293 
Fluctuation of temperature difference 

across PM and housing: 5T d [K Hz" 1 / 2 ] 1.0 X 10~ 5 

Maximum charge build-up: q [C] 10~ 12 

Fluctuation in charge: 8q [C Hz" 1 / 2 ] 6.1 X 10~ 15 

Magnetic susceptibility: \ m 5 X 10~ 5 

Magnetic remanent moment: \M r \ [A m 2 ] 1 X 10~ 7 

Residual gas pressure: P [Pa] 10~ 5 

Electrostatic shielding factors: £ e 10 

Magnetic shielding factors: £ m 1 



Spacecraft (SC) 

Mass: M sc [kg] 350 

Velocity: v [m s" 1 ] 4 X 10 4 

Thrustcr noise [fiN Hz" 1 / 2 ] 10 

Area facing the Sun: A sc [m 2 ] 5 

Fluctuation of temperature in SC: <5T SC [K Hz" 1 / 2 ] 0.2 



Capacitive sensing 

Capacitance: C x [pF] 6 

Capacitance to ground: C g [pF] 6 

Total capacitance: C [pF] (6C 9 ~)36 

Gap: d [mm] 2 

Asymmetry in gap across opposite sides of PM: Ad ]/im] 10 
Average voltage across opposite faces: V XQ (= (V x i + V x 2)/2) [V] 0.1 

Proof mass bias voltage: Vmo [V] 0.6 

Voltage difference between V X(I and voltage to ground: Vo g [V] 0.05 

Voltage difference between opposing faces: Vd [V] 0.01 
Fluctuation of voltage difference across opposite faces: 

SV d [V Hz" 1 / 2 ] 10- 4 

Residual dc bias voltage on electrodes: Vq [V] 0.1 

Loss angle: <5 10 ~ 5 



Magnetic fields 

Local magnetic field: B 3C [T] 8 X 10~ 7 

Local magnetic field gradient: |VB SC | [T m^ 1 ] 3 X 10~ 6 

Fluctuation in local magnetic field: |<5B SC | [T Hz" 1 / 2 ] 1 X 10" 7 

Interplanetary magnetic field: Bi v [T] 1.2 X 10~ 7 

Interplanetary magnetic field gradient: \5Bi p \ [T Hz -1 / 2 ] 4 X 10~ 7 
Gradient of time- varying magnetic field: \V(SB)\ [T m" 1 Hz" 1 / 2 ] 4 X 10" 8 



Cosmic rays 

Impact rate: A [s _1 ] 30 

Mass of positron: m [kg] 1.7 X 10~ 27 

Incident energy: E d [J] (200 MeV =) 3.2 X lO" 11 



Gravity gradients 

Source mass: M dis [kg] 1 in f gg and 0.03 in a gg 

Distance from the source mass: x [m] 0.5 in f gg and 0.05 in o. 

Thermal expansion coefficient: CTE(Aluminium) [K^ 1 ] 2.5 X 10~ 5 



Patch field 

Multiplicative factor: 7 5 
Patch field voltage: V pe [V] 0.1 
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Table 2. Physical constants used in the estimation 



Symbols Values used in the estimations 



Elementary charge 


e 


i.6 xicr 19 c 


Permeability of vacuum 




1.26 xlO" 6 N A -2 


Permittivity of vacuum 




8.9 x 10" 12 F m" 1 


Speed of light in vacuum 


c 


3.0 xlO 8 m s" 1 


Boltzmann constant 


fes 


1.38 xl0~ 23 J K" 1 


Stefan-Boltzmann constant 


(7 


5.7 xl(T 8 W m~ 2 K~ 4 


Gravitational constant 


G 


6.7 X lCT 11 m 3 kg -1 s -2 



Table 3. The estimated values for the proof-mass acceleration disturbances at 
the frequency of 0.1 mHz. PM and rss denote proof mass and root-sum-square, 
respectively. The values are in units of 10 -15 m s — 2 Hz -1 ' 2 . See sections 151 and 
l6l and tables 111 and l2l for notation. 

[10" 15 m s" 2 Hz" 1 / 2 ] Expressions ASTROD I 



(a) PM environmental acceleration disturbances 

/mi Magnetic -^f- | VB SC \8B SC 1.2 

/ma Magnetic y^n.\yB sc \SB ip 3.4 

fm3 Magnetic -^—\M r \\v(SB)\ 1.6 

fm Magnetic ^JsqB ip 0.0017 

/L2 Magnetic ^jJL q SB ip 0.91 

fc Cosmic rays . imE d x 0.0015 



2 \/W (3fcsTpmjv)1/4 a74 



f rg Residual gas 

f re Radiometric pP 0.17 

fog Outgassing 10/ re 1.7 

ftr Thermal radiation f^T7 T p 5T d 013 

fgg Gravity gradients 2G ^ d " CTE ■ 6T SC 2.7 

rss(/ml /ffff) — fnep 5.2 

(b)PM sensor back-action acceleration disturbances 



fbl 




SV d X V g 


^d^V0gSV d 


1.4 


tb2 




8V d X q 




4.8 


fb3 




Sq X V d 




2.9 


fb4 




5q X q 


rrtpd C 2 d * ^ 


0.040 


fba 




Readout electronics 


(see table 1 of 13711 


1.8 


fpe 




Patch fields 




29 


fdl 




Dielectric losses 


1.6 


rss(/f,i - 


fdl ) — fnbp 






30 


rss(/ m i 


- fdl) = fnp 






30 
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Lower susceptibility may be achievable by controlling the manufacturing process of 
the alloy. The requirement for LISA is Xm = 3 x 10 -6 HJ. We use a somewhat 
moderate value of \ rn = 5 x 10 -5 . Assuming £ m = 1 as per ^5]^], we obtain f ml 
= 1.2 x lO-^ms^Hz-Va and f m2 = 3.4 x 10~ 15 ms" 2 Hz" 1 / 2 at O.lmHz. 

Acceleration disturbance due to the magnetic remanent moment is given by f m 3 
listed in table [31 (equation 4(c) of |18|L where A(SB) is time- varying magnetic field 
gradients. According to measurements by Gill et al., magnetic remanent moment of 
three 4.2-cm 70/30 Au-Pt alloy cubes (one of them was 73/27 Au-Pt alloy) ranged from 
3.5 x 10 _9 Am 2 kg _1 to 3.1 x 10~ 8 Am 2 kg _1 [22]. By scaling the largest measured 
remanent moment by weight, it is 5.4 x 10 -8 Am 2 for a 1.75-kg proof mass. By using 
a somewhat relaxed value of 1 x 10~ 7 Am 2 and |V(£B)| = 4 x 10~ 8 T m _1 Hz~ 1/2 
|18| . which is a factor of four higher than the assumed requirement for LISA in |23| . 
we obtain / m 3 = 1.6 x 10~ 15 ms -2 Hz -1 / 2 at O.lmHz. In the error estimates for LISA 
23 , magnetic remanent moment of 2 x 10~ 8 Am 2 is used. 

5. 2. Lorentz force 

The proof mass in orbit would be charged up by cosmic-ray impacts |25l El • As the 
charged proof mass moves through the interplanetary magnetic field (B ip ) with a 
velocity v of about 4 x 10 4 ms _1 , it experiences the Lorentz force: Fl = qv x Bi P , 
where q is the built-up charge. Acceleration disturbances due to the fluctuation of 
the charge build-up (Sq) in the proof mass and of the average interplanetary magnetic 
field are given by fax and fh2 (see table |3J), respectively, where £ e is an electrostatic 
shielding factor; in the rest frame of the proof mass, the motion of the proof mass 
through the interplanetary magnetic field generates an electrostatic field. 

By assuming the Poisson distribution of cosmic-ray impacts, the average charge 
fluctuation spectral density can be defined as Sq(u>) = \J1eqjuj. For the frequency 
[y = (27r) _1 o;) of O.lmHz and the effective charging rate (q) of 288 +es _1 , which 
was estimated for a LISA proof-mass (46-mm cube) by a simulation using GEANT 
4 toolkit ESI, we obtain Sq(uj) = 6.1 x 10~ 15 C Hz^ 1 / 2 . Using this value, £ e = 10 
and B vp = 1.2 x 10" 7 T, we obtain f L1 = 1.7 x 10~ 18 ms~ 2 Hz" 1 / 2 at O.lmHz. The 
maximum value of Bi p at 1 AU is about 3 x 10~ 8 T. We tentatively use four times 
the maximum value as B ip at 0.5 AU, by assuming l/(distance) 2 dependence of B ip $. 
Because the volume of the ASTROD I proof-mass is about 10 % smaller than the LISA 
proof mass, the charging rate for ASTROD I might be smaller than the value for LISA. 
Bao et al. are working on simulations to estimate the charging rates for ASTROD I 
|27| . Stebbins et al. use £ e = 100, which is one order of magnitude larger than the 
value we used here, in the current error estimates for LISA |23j . Taking a nominal 
maximum build-up charge q = 10~ 12 C, which is one order of magnitude larger than 
the value used in the error estimates for LISA |131 1141 12~3*] . we obtain = 9.1 x 
lO-^ms^Hz- 1 ^. 

5.3. Cosmic-ray impacts 

Some cosmic rays get stopped in the proof mass and deposit momentum [2S1 U\- 
Assuming the Poisson distribution of cosmic ray impact, spectral density of momentum 
transfer (p) is 2p 2 X, where A is the fluctuation in the impact rate The acceleration 

tt The solar dipole field at 0.5 AU is less than ~ 1CT 10 T. The main magnetic field at 0.5 AU to 1 AU 
is due to the influence of solar winds which attenuate with (distance) 2 . 
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disturbance due to the fluctuation in the impact rate is given by f c , listed in table 
121 The impact rate was inferred from simulations done for LISA; by adding the 
effects of all stopped particles (protons and helium) and taking into account their 
directions, the acceleration disturbance by momentum transfer was estimated to be « 
2 x 10~ 18 ms -2 Hz -1 / 2 for a LISA proof mass [251 - This corresponds to a disturbance 
due to momentum transfer by protons (mass m =1.7 x 10~ 27 kg), with an impact rate 
of ~ 30s~\ at incident energy E d = 200 MeV (= 3.2 x 10~ n J). Using these values, 
the acceleration disturbance becomes f c — 1.5 x 10~ 18 ms~ 2 Hz -1 / 2 . 

5-4- Residual-gas impacts 

From the kinetic theory, the number of residual-gas molecules (assumed as ideal gas) 
that pass an area (A p ) of the proof mass per second is given by w = nApv/Q, where 



n = P(kgTp) 1 is the number density of the molecules and v = ikpTpm^ is the 

average thermal velocity (P is the pressure of residual gas, kp = 1-38 x 10 -23 JK _1 is 
the Boltzmann constant, Tp is the temperature of the proof-mass housing and mjf = 
4.65 x 10 - 26 kg is the mass of nitrogen molecules). Assuming the Poisson distribution 
of the impact rate, we define the spectral density of fluctuations in w as 5m(uj) = 



Acceleration due to the residual gas impacts is given by 2m^wvm p 1 and 
acceleration disturbance due to fluctuations in the impact rate of residual gas is given 
by f rg listed in tabled! For P = 10~ 5 Pa, A P = 1.75 x 10 -3 m 2 and T P = 293 K, we 
obtain f rg = 7.4 x 10~ 16 ms~ 2 Hz -1 / 2 . Stebbins et al. use 3xlO _6 Pa as the residual 
gas pressure around the proof mass in the error estimates for LISA |23j . 

5.5. Radiometric effect 

Acceleration disturbance due to the radiometric effect (e.g. (2312111201121122) is 
given by f re listed in table [21 where 5Td is fluctuation in temperature difference 
across the proof mass housing. This value has to be estimated by carrying out 
thermal modeling. According to thermal analysis for LISA, temperature fluctuation 
on the optical bench, to which the proof mass housing is mounted, due to power 
dissipation of amplifiers is about 3.0 x 10~ 5 K Hz" 1 / 2 at 1 mHz (table 6.2-28 of |U). 
By assuming that the fluctuation rises as 1/v, the temperature fluctuation of the 
optical bench is ST ob = 3.0 x 10~ 4 K Hz" 1 / 2 at 0.1 mHz. At the frequency of 0.1 mHz 
and higher frequencies, the temperature fluctuation on the optical bench would be 
dominated by the fluctuation in the power dissipation |33| : fluctuation due to solar 
irradiance at 0.1 mHz is 1.1 x lO^KHz" 1 / 2 (table 6.2-16 of 0). The ratio ST ob /ST d 
would range from 30 to 100 (221- By using a value of 30 for the ratio, we obtain 
ST d = 1.0 x lO^KHz- 1 / 2 at 0.1 mHz. By using this value, we obtain f re = 1.7 x 



5.6. Temperature dependent outgassing effect 

Outgassing from walls of the sensor cage is thought to produce greater acceleration 
noise than the radiometric effect |2U • An analysis done for the LISA Technology 
Package (LTP), assuming a simple model of flow circuit with a linear approximation, 
shows that the outgassing effect is nearly 10 times the radiometric effect By 





lO-^ms^Hz- 1 / 2 at 0.1 mHz. 
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using this estimate and the estimate we made for the radiometric effect in the previous 
section, we obtain f og = 1.7 x lfr 15 ms~ 2 Hz" 1 / 2 at O.lmHz for ASTROD I. 

5. 7. Thermal radiation pressure 

By assuming a perfectly reflecting surface of the proof mass, thermal radiation 
pressure produces acceleration disturbance (ft r in table [3J due to fluctuations in the 
temperature difference across the proof-mass housing. In the expression of ft r listed in 
tabled a = 5.7 x 10~ 8 Wm~ 2 K~ 4 is the Stefan-Boltzmann constant, A p is the area 
of the proof mass and T p is the temperature of the proof-mass housing. A factor of 
one-third is multiplied, as done in the estimation for LISA by Schumaker | 13l 114) , as a 
margin for the fact that not all of the radiation momentum is normally incident on the 
proof mass. For the housing temperature of 293 K and the temperature fluctuation 
of 1.0 x lO^KHz- 1 / 2 , we obtain f tr = 1.3 x 1CT 16 ms" 2 Hz" 1 / 2 at O.lmHz. The 
same value was used for the temperature fluctuation in the error estimates for LISA 



5. 8. Gravity gradients due to thermal distortion of the spacecraft 

The temperature fluctuation in solar irradiance would cause fluctuation in distortions 
of the spacecraft: 8xx~ x = CTE ■ |£T SC |, where CTE and 8T SC are the coefficient of 
thermal expansion and temperature fluctuation of the spacecraft, respectively. The 
inherent fluctuation in solar radiation is 5Wq/Wq — AST/T p rs 2.8 x 10~ 3 Iiz -1 / 2 
at O.lmHz ^S]- Therefore, the temperature fluctuation is 0.2 KHz" 1 / 2 for T p = 
293 K. The gravitational disturbance by a 1-kg mass (M^is) separated from the 
proof mass in the sensitive axis by x = 0.5 m of aluminium structure would be 
f gg = 2.7 x 10" 15 ms" 2 Hz- 1 / 2 for ST SC = 0.2 KHz" 1 / 2 and CTE (of Aluminium) 
= 2.5 x 10 _5 K _1 . The same disturbing mass (Mdis = 1kg and x = 0.5 m) is 
assumed in ^2 1141 123| . In reality, the mass to be involved in thermal distortion 
would be much larger, but the influence would be largely canceled because of the 
axial symmetry in the original spacecraft geometry. The inherent fluctuation in solar 
radiation could be reduced largely by thermal shielding. In the estimation for LISA, 
5T SC = 0.004KHz~ 1 / 2 is used For a more accurate estimate, gravity effects 

by thermal and non-thermal distortions of the spacecraft and the payload have to be 
studied by appropriate modeling. 

6. Proof-mass sensor back-action acceleration disturbances 

The total mechanical energy of the capacitive sensing system can be expressed as (e.g. 
equation (A. 3) of |13j): 



where q is a net charge of the proof mass; C is the sum of the capacitances due 
to the applied voltages on the surrounding electrodes i and the potential to ground 
g : C = ^2^Ci, where I = x\,Xi,y\^y%^Zx,z-2,iQ\ V s is the voltage induced on the 
proof mass by the applied voltages on the electrodes: V s = C _1 ^ The first 

term of W is the total energy done on the proof mass by the applied voltages of the 
surrounding electrodes. The second term is the energy acquired on the proof mass by 




(10) 
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the image charge on the surrounding electrodes. The third term is the energy stored 
on the proof mass by the deposit of the free charge on the proof mass. 

The x-component of the force on the proof mass is given by differentiation of 
equation (|10|) with respect to x. For a simplicity, we assume that neither the free 
charge q nor the potentials Vi have appreciable gradients along the x-axis, as per p|] . 
Acceleration disturbances due to fluctuations in the applied voltages and charge can 
be given by fbi - fbi as listed in table [21 A detailed description on the deviations of 
fbi - fbi is given in Appendix A of |13| . The parameter values used in this section 
and the estimated values are listed in tables ^ and [3] respectively. 

To simplify the analysis, several assumptions were made in the process of deriving 
the expressions for the four classes of acceleration disturbances (fbi - fbi) E3 : (1) Q 
are all comparable with each other in magnitude and on the order of C x ~ 6 pF (= 
e^Ad^ 1 , where eq = 8.9 x 10 _12 Fm _1 is the permittivity of vacuum, A is the area 
of each electrode and d = 2mm is the gap); (2) only the capacitances C x \ and C X 2 
have nonzero gradients along the x-axis: C' xl = —C' x2 , and the gradient of the total 
capacitance C ps 6 C x is: C ps (C x Ad)d~ 2 , where Ad is the gap asymmetry in the x- 
direction; (3) the average of the potentials on opposing faces is same for all three axes 
and expressed as V xo (= (V x i + V x2 )/2 = 0.1 V); (4) the magnitude of the fluctuation 
in potential 8Vi is all identical and take the value of the average fluctuation of all the 
potentials for the three axes and the voltage to ground. We express the fluctuation as 
SV X() . It should be noted that we do not consider any cross-talk effects that arise in 
the sensitive axis due to forces applied to the other degrees of freedom. 

6.1. Fluctuations in voltage imbalance and charge 

fbi is associated only with sensing voltages but not the free charge, and its value is 
fbi ps 1.4 x 10- 15 ms- 2 Hz- 1 / 2 , by assuming V 0g = V Xo - V g = 0.05 V. This value is 
larger than the value used for LISA [22| by a factor of five. Also, we use 5Vd = 1.0 
x 10~ 4 VHz -1 / 2 as fluctuation in voltage imbalance Vd (= V x \ — V x2 — 0.01 V). This 
value (5Vd) is one order of magnitude larger than the value used for LISA [T3lll4l I23| . 
Further, we assume that C g /C pa 1/6, as per |13M14| . 

fb2 and fb3 arise from fluctuations in the force due to the interaction between 
the net free charge and applied sensing voltages. fb% is due to fluctuation in voltage 
imbalance and its value is fb2 ~ —4.8 x 10~ 15 ms -2 Hz -1 / 2 , where the net free charge 
on the proof mass is set to the nominal maximum built-up charge, described earlier. 
fb3 is due to fluctuation in charge and its value is fb3 ~ —2.9 x 10~ 15 ms~ 2 Hz -1 / 2 . 
The value of Vd used here is larger than the one used in the error estimates for LISA 
23 by a factor of two. fbA is associated only with the free charge and its value is fbi 
ps 4.0 x 10~ 17 ms~ 2 Hz -1 / 2 . In the estimation, we assume that the gap asymmetry in 
the x-direction is Ad = 10/Ltm, which is larger than |13U14) by one order of magnitude. 

6.2. Readout electronics 

Readout electronics for the capacitive sensing will be similar to ones studied for LISA 
(e.g. [213 E3 Ell)- For the standard resonant inductive-bridge scheme discussed in 
|37|. the sources of disturbances due to readout electronics can be classified into 
two categories: imperfections in the capacitance bridge and the electric noise in 
the detecting circuit. The former includes fluctuations of inductance imbalance 
(AL/L), fluctuations of mutual inductance imbalance (AM/M), bias oscillator 



Acceleration disturbances and requirements for ASTROD I 



13 



relative amplitude noise (AV/Vmo, where Vmo is the 100 kHz bias voltage capacitively 
applied to the proof mass by injection electrodes on insensitive faces (the y and/or 
z faces) of the proof mass; the sensing electrodes are grounded) and bias oscillator 
phase noise (A<fi). The latter includes current noise and thermal noise. 

Using the expressions for these sources of the disturbances given in table 1 of [23 , 
the main contributions to the acceleration disturbance at O.lmHz for ASTROD I are 
from AV/Vmo (~ 1-5 x 10~ 15 ms~ 2 Hz^ 1 / 2 ) and A(f> (w 9.1 x 10~ 16 m s~ 2 Hz" 1 / 2 ), 
and the thermal noise (~ 3.8 x 10 -17 ms~ 2 Hz -1 / 2 ) and the current noise (« 1.3 
x 10~ 17 ms~ 2 Hz" 1 / 2 ). Contributions from fluctuations of AL/L and AM/M are 
insignificant. In the estimation, we use the following parameters based on experimental 
results reported in 37 : AL/L » 10~ 4 , AM/M «6x 10~ 8 , AV/Vmo ~ lO^Hz- 1 /*, 
A<f> w 5.7 x 10 -4 Hz -1 / 2 and the turn ratio of the transformer n = 1 (40 turns, 
inductance of 5mH and a quality factor of 165), and a somewhat relaxed residual 
imbalance of the bridge pdc ~ Ad/d = 5 x 10 -3 . This relaxation has resulted in 
the dominant contributions from the bias oscillator relative amplitude noise and the 
phase noise. In [37], the thermal noise is dominant as they use pdc ~ 10~ 4 in their 
estimation. The rss of these disturbances is fb a ~ 1-8 x 10~ 15 ms -2 Hz -1 / 2 at 0.1 mHz 
and listed in table OH 

6.3. Patch field voltage 

Even when the sensing voltages are not applied on the electrodes, differences in 
local surface properties of the electrodes and the proof mass could lead a potential 
difference, patch- field voltage, between them j^Hj- The charge fluctuations 8q result 
in acceleration disturbance f pe (table 03) through the patch field [2B]- By taking 
the average patch-filed voltage difference between opposing electrodes as V pe =0.1V 
[E3HI1EB], we obtain f pe = 2.9 x 10~ 14 ms~ 2 Hz" 1 / 2 at 0.1 mHz. This is the dominant 
contribution to the total acceleration disturbance of the ASTROD I proof mass. 
The LISA study team is investigating the possibility of measuring and compensating 
voltage differences across capacitors, during the mission commission process, to 
considerably better than 0.01 V |23| . 

6.4- Dielectric losses 

Dielectric losses are thought to stem from surface contamination of electrodes and 
produces thermal voltage noise 001 EH] : 



Svdiel = \ Uk B T p — — (11) 

V ULj x 

where S is loss angle. The upper limit of 6 is reported to be 10 -5 for Al electrodes 
|41| . For S — 10~ 5 , this voltage noise is about 6.6 fj,Y Hz -1 / 2 at 0.1 mHz and produces 
acceleration disturbance fdi (see table 01 for the expression), through residual dc bias 
voltage on electrodes, in the sensitive axis |40l I38j . By making the same assumption 
of the average potential difference Vo = 0.1 V between a given electrode and the proof 
mass as |3H], we obtain fdi ~ 1.6 x 10~ 15 ms~ 2 Hz" 1 / 2 . In the error estimates for 
LISA [231, Vb = 0.01 V is used. 
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6.5. Summary of the direct acceleration disturbances of the proof mass 

By adding in quadrature, the rss of the proof-mass environmental acceleration 
disturbances (fnep) and the sensor back-action acceleration disturbances (fnbp) are 
5.2 x 10" 15 ms- 2 Hz" 1 / 2 and 3.0 x 1CT 14 ms~ 2 Hz" 1 / 2 , respectively. Therefore, the 
total direct proof-mass acceleration disturbance (fnp) is 3.0 x 10~ 14 ms~ 2 Hz" 1 / 2 at 
0.1 mHz, dominated by the sensor back-action acceleration disturbances. 

7. Proof mass-spacecraft coupling 

The stiffness K is considered to stem from gravity gradients (K gg ), fluctuations in 
sensing capacitance and capacitance gradients (K s i, K s2 and K s s), bias voltage (K S 4), 
patch field voltage (K s5 ) and magnetic field gradients {K m \ and K m2 ). Table^gives 
a summary of estimated values for these sources of stiffness. The expressions used in 
the estimations are briefly described below. A detailed description on deviations of 
the expressions for K gg , K s i, K s2 , and Ks$ is given in Appendix A of |13|tt. 

7.1. Gravity gradients 

For a given disturbing point mass (Mdi S ) at a distance x from the center of mass of 
the proof-mass on the sensitive axis, the amplitude of the acceleration disturbance of 
the proof-mass caused by a positional fluctuation X ps is a gg = K gg X ps , where K gg is 
given m table H 

Making the same assumption of Mdis — 0.03kg and x = 0.05 m as [22], we obtain 
K gg as 3.2 x 10~ 8 s~ 2 . For a more detailed analysis, the identification of the disturbing 
mass is necessary. Gravitational modeling for ASTROD I is in progress |42| . 

This disturbance arises from any positional fluctuation and is different from the 
gravity gradient caused by thermal distortion or motion (f gg ), which was discussed 
earlier. 

7.2. Fluctuations in capacitive sensing 

The expressions for K s i, K s i and K s3 can be obtained in the similar way as fbi, fb2, 
fb3 and fb4 in the previous section, under the following assumptions fluctuations 
in the capacitances C x \ and C x2 and their derivatives produce disturbing forces in 
the x-direction, but not the fluctuations in other capacitances or their derivatives; we 
ignore the cross-coupling effects. Also, we assume 8C X \ = —5C X 2 ~ (C x /d)Sx and 
SC' xl = SC' x2 « (C x /d 2 )6x. 

K s \ is due to the fluctuations in the Coulomb interaction between the charged 
proof mass and the image charges on the surrounding electrodes. K s2 arises from 
interaction between the net free charge q on the proof-mass and the average electrode 
voltages. K s3 is due to the applied voltages across electrodes and the voltage difference 
across opposite electrodes. The estimated values for K s i, K s2 and K s ^ are given in 
table H 

7.3. Bias voltage 

By assuming that the sensing electrodes are grounded, the dominant term of the 
readout stiffness along the sensitive axis is given by K S 4 |2Hj, listed in table 0] 

tt-^«5 i n text is noted as K s 4 in ll.'il . 
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Table 4. The proof-mass stiffness in units of 10 9 s 2 (see sections |§] and and 
tables [T1 and |21 for notation.) 

[10 -9 s" 2 ] Expressions ASTROD I 

K gg Gravity gradients 2 GI ^$ ZB 32 

K sl Image charges - ^ (^) q 2 0.66 

K s2 q xV 0g (LiJ (%) (^) 9 y 09 0.40 



v; 2 9 !> 0.095 



^3 Applied voltages ^ j + i) V 2 + (^) 

K s4 Bias voltage -^3? V M0 3.1 X 10 2 

K a5 Patch fields 7 (-^ \ p % 1.2 

ifml Induced magnetic moments 22-™. { I VB SC | 2 + B sc \ V 2 B SC \ ) 9.0 X 10~ 5 

K m 2 Magnetic remanent moments — |M r .||V 2 -B sc | 6.9 X 10~ 4 

rss(i<" 99 - K m2 ) = K 3.1 x 10 2 



Assuming, as per that the proof mass is biased to Vmo — 0.6V, K s a ~ 3.1 

x 10" 7 s" 2 . 

7.4- Patch field voltage 

By assuming a nominal patch-field voltage of V pe = 0.1 V |13l 1141 155] and an overall 
multiplicative factor of 7 = 5 as per | 131 IT4*| , the contribution due to the patch field 
to the proof mass-spacecraft coupling is K s5 w 1.2 x 10 _9 s -2 (table 

7. 5. Magnetic field gradients 

The magnetic stiffness is given by K m = —V[V(M P • B)\ ^8]. The expressions 
for dominant terms of the stiffness due to the induced magnetic moment and the 
remanent moment are given by K ml and K m2 (as listed in tabled}, respectively, 
where |V 2 B SC | w 12B sc r m 2 = 1.7 x 10~ 5 Tm~ 2 for r m = 0.75m. Their contributions 
to the total stiffness are insignificant. 

7. 6. Summary of the estimated values of stiffness 

These contributions to the coupling constant K are summarized in table 21 The rss 
of the coupling constant is 3.1 x 10~ 7 s~ 2 . This is slightly below the requirement for 
the total stiffness in LISA (4 x 10~ 7 s~ 2 |55]1. 

8. Requirements for the readout sensitivity and spacecraft control-loop 
gain 

We have estimated values for the coupling constant K, the direct spacecraft 
acceleration disturbance f ns and the direct proof-mass acceleration disturbance f np . 
By using the expression for the total acceleration disturbance of f p (equation (JSJ)), 
we infer the requirements for the readout sensitivity X nr and the spacecraft control- 
loop gain u. In this process, we allocate an identical magnitude f a to each term 
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of the expression; ft = /£/3. For ASTROD I, the noise goal is 10~ 13 ms~ 2 Hz" 1 / 2 
at O.lmHz and, therefore, f a = 5.8 x 10 -14 ms" 2 Hz -1 / 2 . The second term of the 
expression is f np , which was estimated to be 3.0 x 10~ 14 ms~ 2 Hz~ 4 / 2 (table 0; it is 
about a factor of 2 smaller than the allocated requirement. With the estimated total 
stiffness K = 3.1 X 10~ 7 s~ 2 , we obtain X nr < 1.9 x 10 _7 mHz" 1 / 2 from the first term 
and u > 3.8 x 10 5 from the last term of the expression. 

From figure one can see that the acceleration noise spectral density 
requirements for ASTROD I take its lowest value of about 0.4 x 10~ 13 ms~ 2 Hz -1 / 2 
at 0.3 mHz. At this frequency, f a becomes 2.3 x 10~ 14 ms~ 2 Hz -1 / 2 . Therefore, the 
requirement for X nr becomes more stringent at 0.3 mHz: X nr < 7.4 x lO^mHz -1 / 2 . 
As for the second term of f p , f np is smaller at higher frequencies ^J^l- Our estimate 
of fnp is dominated by the contribution from f pe , which scales as v . Therefore, at 
0.3 mHz, f np would be ~ 1 x 10 _14 ms _2 Hz, which is smaller than f a at 0.3 mHz. 
The last term scales as w~ 2 , and f ns is expected to be smaller at higher frequencies 
because of the v~ x l z dependence of the fractional fluctuation in solar irradiance (figure 
6 of |15p. Therefore, by assuming the same level of the thruster noise, the requirement 
for u at 0.3 mHz would be less stringent than that at O.lmHz by a factor of ~ 4. 

In summary, the requirements for the readout sensitivity and the control loop 
gain for ASTROD I are X nr < 7.4 x 10~ 8 mHz- 1 / 2 and u > 3.8 x 10 5 , respectively. 

9. Comparison with LISA 

Main relaxed parameter- values are listed in table El in comparison with LISA. The 
values for LISA are quoted from the current error estimates by Stebbins et al. |23| . 
except the thruster noise quoted from Recently the LISA requirements for the 
thruster noise and the residual gas pressure have been relaxed further from the values 
given in table (this is a comment from one of the referees). 



Table 5. Relaxed parameter values in comparison with LISA 


v = 0.1 mHz 


ASTROD I 


LISA 


Maximum charge build-up: q [C] 


lO" 12 


lO" 13 


Magnetic susceptibility: \ m 


5 x 10~ 5 


3 x 10~ 6 


Magnetic remanent moment: \M r \ [A m 2 ] 


1 x 10~ 7 


2 x 10" 8 


Residual gas pressure: P [Pa] 


10" 5 


3x 10" 6b 


Electrostatic shielding factors: £ e 


10 


100 


Thruster noise \pN Hz -1 / 2 ] 


10 


a few 


Fluctuation of temperature in SC: ST SC [K Hz -1 ' 2 ] 


0.2 


0.004 


Voltage difference between average voltage across 






opposite faces and voltage to ground: Vo g [V] 


0.05 


0.01 


Voltage difference between opposing faces: Vd [V] 


0.01 


0.005 


Fluctuation of voltage difference across opposite faces: 






SV d [V Hz -1 / 2 ] 


10" 4 


10~ 5 


Residual dc bias voltage on electrodes: Vq [V] 


0.1 


0.01 



" According to one of the referees, this LISA requirement has been relaxed to 
10" 5 Pa. 
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Table 6. Acceleration disturbances and requirements. PM and SC denote proof- 
mass and spacecraft, respectively. 

ASTROD I 

(a) Estimated total acceleration disturbance (at 0.1 mHz): 
f p [m s" 2 Hz-V2] s.7 x 10" 14 

(b) Estimated contributions to f p (at 0.1 mHz): 

f p ~ X nr (-K) + f np + (f ns + TN^Kuj^u- 1 
Direct PM acceleration disturbances: 
f„p [m s- 2 Hz-V2] 3.0 x 10~ 14 

Direct SC acceleration disturbances: 
Us + TN t [m s- 2 Hz-V 2 ] 2.8 x 10~ 8 

PM-SC spring constant: 

K [s" 2 ] 3.1 x 10~ 7 

(c) Inferred requirements 

Position readout noise: X nr [m Hz" 1 / 2 ] 7.4 X 10~ 8 
Control-loop gain: u 3.8 X 10 5 



10. Summary, discussion and conclusions 

We have estimated the spacecraft acceleration disturbance f ns (section 0J, the proof- 
mass acceleration disturbances f np (sections [5] and EJ and the stiffness K between 
the spacecraft and the proof mass (section 0) . By using the expression (fHJ , we have 
inferred the requirements for the readout sensitivity X nr and the control-loop gain u 
(section - 

Table provides a summary of the estimated acceleration disturbances and 
the stiffness (section (b)), and the requirements (section (c)). The estimated total 
acceleration disturbance f p at 0.1 mHz (section (a) of table is about 13 % less than 
the noise goal of 1CP 13 ms~ 2 Hz~ 1//2 . We have compared the parameter values used 
in the estimation with LISA in sectional 

The total direct acceleration disturbance of the proof mass (f n p) at 0.1 mHz was 
estimated to be nearly a factor of two smaller than the requirement. This ^50% 
margin may be allocated for unknown disturbances or disturbances that would arise 
but have not been studied yet. These unestimated disturbances would be originated 
from, for instance, cross-talks in the capacitive sensing and magnetic damping of 
the proof mass. An estimate of acceleration disturbance due to magnetic damping 
for LISA is about 2 x 10~ 16 ms~ 2 Hz" 1 / 2 at 0.1 mHz 23 . The contribution from the 
magnetic damping effect to ASTROD I would be in the similar order and insignificant. 

The sensor back-action acceleration disturbances can be reduced by increasing the 
magnitude of the gap d. The total stiffness K would be also reduced, for example, by 
a factor of four by changing the gap to 4 mm. The optimum design for the capacitive 
sensing is to be discussed based on results from the ongoing laboratory torsion balance 
experiment for ASTROD I gSj. 

Parameter values we used in this paper are mainly based on the results of studies 
done for LISA and LISA Pathfinder. This may be sufficient for the preliminary 
estimation. More accurate estimation would be obtained by carrying out the following 
works dedicated for ASTROD I: (a) modeling local magnetic fields of the spacecraft, 
(b) estimating the effective charging rate of the proof-mass, (c) estimating the cosmic- 
ray impact rate of the proof-mass, (d) estimating the micrometeorite impact effects, 
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(e) thermal modeling of the proof- mass housing and the spacecraft, (f) gravitational 
modeling that includes thermal and non-thermal deformation of the spacecraft and 
the payload, (g) electrostatic modeling for the capacitive sensors and (h) estimating 
environmental factors (such as the interplanetary magnetic field, solar wind, solar 
radiation and cosmic rays) in the varying orbit (0.5 AU to 1AU). Simulations to 
estimate charging rates for ASTROD I are in progress |27| . 

We have tentatively estimated acceleration disturbances for ASTROD I. This 
work has allowed us to set preliminary requirements for ASTROD I. To improve 
the current estimation, the disturbances that have not been studied yet have to 
be included and more detailed modeling works are necessary for ASTROD I. In 
comparison with LISA, requirements for ASTROD I can be largely relaxed. This 
will make the technological developments for ASTROD I less demanding to meet the 
drag- free requirements. 
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